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ABSTRACT 

Context. Betelgeuse is one the largest stars in the sky in terms of angular diameter Structures on the stellar photosphere have been 
detected in the visible and near-infrared as well as a compact molecular environment called the MOLsphere. Mid-infrared observations 
Ph \ have revealed the nature of some of the molecules in the MOLsphere, some being the precursor of dust. 

Aims. Betelgeuse is an excellent candidate to understand the process of mass loss in red supergiants. Using diffraction-limited adaptive 
optics (AO) in the near-infrared, we probe the photosphere and close environment of Betelgeuse to study the wavelength dependence 
of its extension, and to search for asymmetries. 

Methods. We obtained AO images with the VLT/NACO instrument, taking advantage of the "cube" mode of the CONICA camera 
to record separately a large number of short-exposure frames. This allowed us to adopt a "lucky imaging" approach for the data 
1 reduction, and obtain diffraction-limited images over the spectral range 1.04-2.17 yt/m in 10 narrow-band filters. 

Results. In all filters, the photosphere of Betelgeuse appears partly resolved. We identify an asymmetric envelope around the star, 
with in particular a relatively bright "plume" extending in the southwestern quadrant up to a radius of approximately six times the 
photosphere. The CN molecule provides an excellent match to the 1.09 iim bandhead in absorption in front of the stellar photosphere, 
but the emission spectrum of the plume is more difficult to interpret. 

Conclusions. Our AO images show that the envelope surrounding Betelgeuse has a complex and irregular structure. We propose that 
the southwestern plume is linked either to the presence of a convective hot spot on the photosphere, or to the rotation of the star. 

co : 

i Key words. Stars: individual: Betelgeuse; Stars: imaging; Stars: supergiants; Stars: circumstellar matter; Techniques: high angular 

qq . resolution; Methods: observational 

■ 1. Introduction needs to be invoked. A careful scrutiny of the object from the 
' photospheric scale up to the distance in its circumstellar enve- 
; Betelgeuse (a Ori, HD 39801, HR 2061) is one of the brightest lope (CSE) where dust can conde nse is the path to understanding 
. . . stars at infrared wavelengths and has one of the largest apparent this prototypical and enigmatic star. It is also known from inter- 
photosphenc diameters (-44 mas in the K band, from Perrin et f er ometric observations in the thermal infrared domain by Bester 
^ al.|2004D. It is a red supergiant of spectral type M2Iab with irreg- et al ( 1996) that Betelgeuse experiences episodic mass loss. Our 
5_J ular flux variations. Because of its brightness and of its angular new NAC0 observations aim at probing the close environment 
Ctf ; size, it has been an easy target for interferometers. It can even be of Betelgeuse to etablish a link between the photosphere and 
resolved with large telescopes at short wavelengths, as demon- mass loss After a presentation of our observations in Sect. [2 we 
strated by Gilliland & Dupree U996}, who resolved its extended detail our ana i ysis of the NACO images in Sect. [3 and discuss 
chromosphere using the Hubble Space Telescope at ultraviolet observed features in Sect H 
wavelengths. Several mysteries remain to be solved for such an 
object among which the structure of its convection and the mech- 
anism of its mass loss. Similarities with Mira stars have been 2. Observations and data reduction 
noted, including a significant thermal infrared excess, indicative 21 lnstrumental setup and observation log 
of the presence or circumstellar dust. An important difference is 

that Betelgeuse is not subject to regular large amplitude pulsa- We observed Betelgeuse on the nights of 2 and 3 January 2009 

tions which could levitate material high enough above the atmo- using the Nasmyth Adaptive Optics System (NAOS, Rousset 

sphere where dust could form and be blown away by radiation et al. 2003) of the Very Large Telescope (VLT), coupled to 

pressure. The reason for the extension of the atmosphere is still the CONICA infrared camera (Lenzen et al. 1998 ), abbrevi- 

unclear and another mechanism for dust formation and mass loss ated as NACO. NAOS is equipped with a tip-tilt mirror and a 

deformable mirror controlled by 185 actuators, as well as two 



Send offprint requests to: P. Kervella wavefront sensors: one for visible light and one for the infrared 

* Based on observations made with ESO telescopes at Paranal domain. For our observations, we exclusively used the visible 

Observatory, under ESO program 082.D-0172(A). light wavefront sensor. The detector is a 1024 X 1024 pixel 
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Fig. 1. Transmissions of the NACO narrow-band filters (black 
curves) and of the atmosphere (from Lord 1992, thin grey curve). 

at the Nasmyth focus of the Unit Telescope 4 (Yepun), located 
on the eastern side of the VLT platform. The CONICA camera 
was read using the newly implemented cube mode that allows to 
record very short exposures (down to 7.2 ms) over part of the de- 
tector. As discussed in Sect. 12.21 this new capability is essential 
to reach truly diffraction limited imaging. 

We selected the smallest available pixel scale of 13.26 ± 
0.03mas/pix (Masciadri et al. 120031 Table [Q gives the list of 
the observations of Betelgeuse and the PSF calibrator stars 
Aldebaran (HD 29139), 310ri (HD 36167) and £Phe (HD 
9362). Due to the brightness of Betelgeuse and Aldebaran, we 
employed narrow- and intermediate-band filters together with a 
neutral density filter (labeled "ND2 .short") with a transmission 
of about 1.5%. The observations of our fainter, secondary PSF 
calibrators (3 1 Ori and 6 Phe) were obtained without neutral den- 
sity filter. The transmission curves of the 10 narrow-band filters 
we used for our observations are shown in Fig.Q] They are also 
available through the NACO instrument web page at ESC0. 

2.2. Data processing 

The individual raw images were pre-processed (bias subtraction, 
flat-fielding, bad pixel masking) using the YorickB and IRAlQ 
software packages in a standard way, except that we did not sub- 
tract the negligible sky background. 

The quality of adaptive optics images is affected by a resid- 
ual jitter of the star on the detector ("tip-tilt"), and by high-order 
residual wavefront distorsions. The star jitter is most probably 
caused by a combination of instrumental vibrations and atmo- 
spheric tip-tilt residuals, while the high-order distorsions are 
simply due to the limitations of the instrument performance. In 
order to obtain truly diffraction-limited image quality, especially 
at short wavelengths in the / band, it is essential to eliminate 
the perturbations from these two effects. Taking advantage of 
the newly implemented cube mode of CONICA, we obtained 
extremely short exposures that allowed us to freeze the image 

1 http://www.eso.org/sci/facilities/paranal/instruments/naco/ 

2 http://yorick.sourceforge.net/ 

3 IRAF is distributed by the NOAO, which are operated by the 
Association of Universities for Research in Astronomy, Inc., under co- 
operative agreement with the National Science Foundation. 



jitter. After a precentering at the integer pixel level, the images 
were sorted based on their maximum intensity, used as a proxy 
of the Strehl ratio. The 10% best images of each cube were then 
resampled up by a factor 4 using a cubic spline interpolation, 
and the star image was subsequently centered using gaussian fit- 
ting. The resulting cubes were eventually averaged to obtain the 
master images of each star used in the image analysis process 
described below. The initial oversampling of the image followed 
by the precise centering of the images gives an improved sam- 
pling of the image, as the residual atmospheric "drizzle" allows 
to retrieve sub-pixel sampling. This is particularly important in 
the / band, where the resolution of the telescope corresponds to 
only two NACO pixels. The selection of the best frames allowed 
us to reject the frames affected by the residual high-order dis- 
torsions. The result of this processing is presented in Fig. [2] for 
Betelgeuse and its primary PSF calibrator Aldebaran. 

2.3. Efficiency of NACO full-pupil cube vs. standard mode 

In order to check the efficiency of the image recentering and 
selection, we processed a series of 10 data cubes obtained on 
the PSF calibrator Aldebaran (#001 to #010 in Table© with and 
without recentering the individual frames. The resulting images 
for three wavelengths are shown in Fig. [3] Already clearly visible 
on these images, we hereafter quantify the gain in image quality 
by processing the same data cubes using three methods: 

- recentering, selection of the best 10% frames based on the 
peak flux in the image and co-addition, 

- recentering and co-addition of all frames (no selection), 

- coaddition without recentering of all frames (equivalent to a 
single long exposure of 72 s). 

From this exercise, it appears that the recentering of the 
frames before co-addition results in a spectacular increase of the 
Strehl ratio and a very significant decrease of the full width at 
half maximum (FWHM) of the star image (Fig.©. As expected, 
this effect is stronger at shorter wavelengths, where the adaptive 
optics performance is comparatively lower. In the 1.04-1. 09 yum 
range, the improvement reaches a factor three in Strehl ratio and 
a factor two in FWHM. The improvement in FWHM brought 
by the selection of the frames with the best Strehl appears very 
small compared to the effect of the recentering. Selecting the 
10% best frames improves the average Strehl by approximately 
10% (relatively to centering without selection) for the shortest 
wavelengths. For our very bright targets, we chose to keep this 
selection criterion as the signal-to-noise ratio of the average im- 
ages is essentially limited by the AO residuals, and not by the 
number of frames. 

As shown in Fig. |H we can also compare the image FWHM 
with the theoretical diffraction limit of the telescope. Using a 
very simplified model of an 8 m primary mirror with a cen- 
tral 1.12 m obscuration, and a 20 mas uniform disk (UD) for 
Aldebaran, we obtained the theoretical FWHM shown as a thick 
grey curve. We chose this UD angular size as the angular di- 
ameter of Aldebaran was measured in the K band by Richichi & 
Roccatagliata (120051 who found 6u D (K) = 19.96+0.03 mas and 
LD (K) = 20.58 ± 0.03 mas. As a remark, the FWHM of our im- 
ages is measured over a large sub-window that includes part of 
the adaptive optics photometric incoherent residuals (also called 
"pedestal"). The resulting values are thus slightly overestimated, 
but have the advantage of being directly comparable between 
real and theoretical images. 

From these results, the frame-by-frame recentering of adap- 
tive optics observations appears as a very efficient way to im- 
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Fig. 3. Average frames of the NACO cubes #001 (1 .04 /mi), #007 
(1.64 /mi) and #010 (2.11 jam) obtained on the PSF calibrator 
Aldebaran without centering or frame selection (left) and with 
centering and 10% best frame selection (right). The white disk 
in the 1 .04 /mi images represents the angular size of Aldebaran 's 
photosphere (20 mas). The color scales are proportional to the 
square root of the flux, and are identical for the two images at 
each wavelength. 



prove the angular resolution compared to single long exposures. 
The spectacular increase in Strehl ratio observed at short wave- 
lengths should also result in an improved sensitivity for faint 
object, but the necessity of combining a large number of frames 
results in a proportionally higher detector readout noise. For this 
reason, this data reduction strategy may not be optimal for K 
band observations where the gain is "only" 50% in Strehl ra- 
tio. However, in the / band, whenever a bright enough source is 
available to recenter the frames, the cube mode of NACO with 
short exposures will in most cases be more sensitive than classi- 
cal long exposures. 



3. Image analysis 

3.1. Point spread function 

We obtained images of two unresolved and relatively fainter 
stars (3 1 Ori and 6 Phe, #022-025 and #05 1-054 in TableQ} with- 
out neutral density filter, to check the influence of the ND fil- 
ter on the PSF shape. We do not find evidence of an effect of 
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Fig. 4. Upper panel: Relative Strehl ratio of the recentered and 
co-added images compared to the classic mode of NACO, with- 
out frame selection (dashed curve) and keeping only the 10% 
best peak flux frames (thin solid curve). Lower panel: Full width 
at half maximum of the NACO images of Aldebaran for the 
classic mode (black curve), frame centering without selection 
(dashed curve) and 10% best frames selection and centering 
(thin solid curve), compared to the theoretical diffraction pattern 
for a UT observing a uniform disk of 20 mas (grey curve). 



this filter on the PSF shape. However, considering the faint- 
ness of the CSE of Betelgeuse, we chose to strictly keep the 
same instrumental setup between Betelgeuse and the PSF to 
avoid any contamination in the deconvolved images. We there- 
fore used Aldebaran as the PSF reference. To retrieve the true 
PSF of the telescope from the images of Aldebaran, we took 
into account the angular extension of the star. Although this 
star is significantly smaller (Old = 20.6 mas) than Betelgeuse 
(#ld = 43.7 mas), it is not negligible compared to the reso- 
lution of the telescope in the / band. We therefore computed 
synthetic images of Aldebaran using its known LD size and 
model limb darkening coefficients in the J, H, and K bands from 
Claret's (12000b tables using Teff = 4 000K and log g = 1.5. We 
then computed a classical Lucy-Richardson (L-R) deconvolution 
of the NACO images of Aldebaran (single iteration) using these 
synthetic images as the "dirty beams". This procedure gives us 
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Fig. 5. Calibrated photometry of Betelgeuse. The dots represent 
our photometry in the 10 narrow-band filters of NACO. The thin 
black curve represents the model spectrum of a T Q s = 3 500 K 
normal supergiant (Castelli & Kurucz 2003 ) normalized to a 
limb darkened disk angular diameter of 43.7 mas. The grey curve 
is a spectrum of Betelgeuse from the IRTF database (Cushing et 
al. l2U05l) . scaled by a factor 1.2. 



the true PSF of the telescope for our ten narrow-band filters, i.e. 
the image of a point source. 

3.2. Photometry of Betelgeuse 

For the photometric calibration of our images, we used a clas- 
sical aperture photometry approach using Aldebaran as the cali- 
brator. As we used non-standard narrow-band filters, we have to 
refer our measurements to a synthetic spectral energy distribu- 
tion (SED) of Aldebaran, normalized by its interferometrically 
measured angular diameter. 

The shape of the SED of Aldebaran was read from the the- 
oretical spectra table by Castelli & Kurucz (2003 ), for an effec- 
tive temperature of 4 000 K, solar metallicity, W = 2 km.s" 1 and 
logg =1.5 (parameters approximated from Cayrel de Strobel et 
al. 120011) . We then used the limb darkened disk angular diam- 
eter of 20.6 mas measured by Richichi & Roccatagliata (2005 ) 
to normalize this synthetic spectrum. We checked that the ex- 
isting broadband photometric measurements are consistent with 
the resulting SED, but we did not used them to fit the SED. 
Considering the true profiles of the NACO filters we used for 
our observations (Fig.[T]), we integrated the flux from Aldebaran 
in W.m" 2 .yum _1 in each band. We then normalized the flux in- 
tegrated on the Betelgeuse images by the flux measured on 
Aldebaran using the same relatively narrow aperture (50 to 
70 mas in radius, slightly increasing with wavelength). This ap- 
proach allows us to integrate the full flux from the resolved 
Betelgeuse disk, while avoiding the adaptive optics halo. We did 
not consider the airmass difference between the two stars and the 
negligible infrared background. 

The resulting photometric measurements are shown in Fig. [5] 
for the ten narrow-band flters. In this figure, we plotted the syn- 
thetic spectrum of a normal supergiant with T G s = 3 500 K, solar 
metallicity, Vt = 2 km.s -1 and logg = 0.0 from the Castelli & 
Kurucz (12003b library. We normalized the expected flux using 
the 43.7 mas angular diameter measured by Perrin et al. (12004b . 



The overall agreement with our photometry is satisfactory, but a 
strong absorption feature around 1.09-1. 24 jjm is present in our 
data and absent from the model spectrum. We also plo tted in 
Fig. [5] a calibrated spectrum of Betelgeuse (Cushing et al. 12005b 
taken from the IRTF databasfl scaled by a factor 1.2 to match 
the brightness of Betelgeuse at the epoch of our observations. 
This spectrum also shows the absorption feature that is particu- 
larly strong at 1.09/zm. 

3.3. Image deconvolution 

We deconvolved the images of Betelgeuse using the PSF images 
(Sect. [XT} as the dirty beams and the Lucy-Richardson (L-R) al- 
gorithm. The images corresponding to the three imaging epochs 
of Betelgeuse were processed separately with the relevant PSF 
calibrator observation. The resulting average deconvolved im- 
ages in the ten NACO narrow-band filters are presented in Fig.O 
We stopped the L-R deconvolution after 30 iterations, as the re- 
sulting images in the K band (2.12 and 2.17/mi) show a full 
width at half maximum of the photosphere of Betelgeuse com- 
patible with the LD size measured by interferometry (Fig. [51 left 
column), while no significant deconvolution artefact is present. 

In this process, we took into account the difference in ro- 
tational position of the pupil of the telescope between the two 
stars. We rotated the images of Aldebaran obtained in each filter 
to precisely match the telescope pupil position angle (hereafter 
abridged as "PA") of the corresponding images of Betelgeuse. 
In other words, the secondary mirror support "spikes" were 
brought to the same PA in Aldebaran 's and Betelgeuse 's im- 
ages. This allowed us to avoid artefacts in the deconvolution 
of the Betelgeuse images. Although the presence of the spikes 
has a noticeable effect on the deconvolved images, we checked 
that the straight deconvolution of the Betelgeuse images by the 
Aldebaran images (without rotation) gives very similar results. 
Photometry on deconvolved images is a tricky task, but McNeil 
& Moody (2005 ) showed that the Lucy-Richardson algorithm 
preserves the photometric accuracy relatively well compared to 
other classical algorithms. We therefore normalized photometri- 
cally the deconvolved images presented in Fig. 0(407 x 407 mas 
field of view) by the total flux of Betelgeuse in each filter (in 
W.m" 2 .yum _1 ) as estimated in Sect. 13.21 and by the size of the re- 
sampled image pixel (3.315 mas/pix). The resulting physical unit 
of the images is therefore W.m" 2 .yum _1 .sr" 1 (surface brightness). 

3.4. Circumstellar envelope shape 

In the deconvolved images presented in Fig. [6l we observe the 
photosphere of Betelgeuse as a disk whose size appears mostly 
constant with wavelength (left columns). We also detect an ex- 
tension of its envelope in the southwestern quadrant, at a PA 
of 200 ± 20 degrees (counted positively east of North), that we 
will subsequently refer to as a "plume". The bright part of this 
plume extends up to a minimum radius of six times the photo- 
sphere of the star 90 - 100 mas), and is already visible in the 
non-deconvolved images of Betelgeuse. At a distance of 4 R*, its 
surface brightness in the 1.09 yum band is « 0.4% of the center of 
the stellar disk, and falls down to « 0.1% at 6 R* (Fig. [8]). Two 
other faint extensions appear in the northeastern (PA ^75°) and 
northwestern (PA « 310°) quadrants. Although faint, their pres- 
ence at all ten wavelengths (processed independently) indicates 
that they are likely real features. 



4 http://irtfweb.ifa.hawaii.edu/~spex/IRTF_Spectral_Library/ 
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Fig. 7. Left panel: Positions of the 3D pseudo-slit, plotted over 
the image of Betelgeuse at 1.04yum. The color scale is logarith- 
mic to better show the extensions of the CSE. The central circle 
represents the photosphere of the star. The graduations from 
to 6 are distances from the photosphere center in stellar radii for 
the spectra shown in Fig.[U Right panel: Reconstructed long-slit 
spectrum of Betelgeuse (linear interpolation from the pseudo-slit 
cut in the Betelgeuse image cube). The color scale is logarith- 
mic. The vertical lines indicate the central wavelengths of the 
NACO narrow-band filters (see also Fig. [T]), and the 43.7 mas 
photospheric angular diameter is shown with green lines. 



3.5. Photospheric and CSE spectral energy distributions 

Based on the deconvolved and photometrically calibrated im- 
ages, we can analyze locally the SED of the CSE of Betelgeuse. 
Fig.|7]shows the pseudo-slit we used for this analysis, positioned 
in a north- south direction to include the brightest plume. The 
right part of this Figure shows the average surface brightness 
spectrum of Betelgeuse and its CSE over the pseudo-slit. The 
individual spectra observed at different distances from the cen- 
ter of the disk of Betelgeuse along the pseudo-slit are shown in 
Fig. [8] (they correspond to horizontal cuts in the right part of 
Fig. [7]). We used the RMS dispersion of the pseudo-slit measure- 
ments on the three deconvolved and photometrically calibrated 
images of Betelgeuse obtained in each filter to estimate the sys- 
tematic uncertainty of our measurements. This gave us the er- 
ror bars shown in Fig. [8] The two fainter plumes visible in the 
Betelgeuse images in the northeastern and northwestern quad- 
rants have similar spectra as the southern plume, with a maxi- 
mum flux in the / band between 1.08 and 1.24 /mi. 



4. Discussion 

4.1. Identification of the spectral features 

Both the photospheric spectrum and the spectra of the plume at 
different radii show a spectral feature from roughly 1 . 1 jim up to 
1.3yUirQ. Against the photosphere, it is seen in absorption, and 
without a strong background source, i.e. offset from the central 
star but at the location of the plume, it is observed in emission. 

The near-IR spectral range of red supergiant stars such as 
a Ori and \i Cep is the topic of a long-standing debate: several 
features have alternately been attributed to CN or H 2 (Woolf et 
al. [T964l Wing & Spinrad IT9701 Tsuii [2000b . The current con- 
sensus is that features due to photospheric CN are present, as 



5 As there is a gap in the NACO spectral coverage from 1.16 to 
1.22yum (see Fig.[T}, it is conceivable that we are in fact dealing with 
two separate features. 
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Fig. 8. Surface brightness spectra of the close environment of 
Betelgeuse, averaged over the pseudo-slit width (44 mas). The 
vertical scales are in W.m~ 2 ./mr 1 .sr~ 1 , and the horizontal scale 
is the wavelength in yum. The radius (in R*) from the center of the 
photosphere is indicated in the upper right corner of each plot. 



predicted by the self-consistent hydrostatic models, and that wa- 
ter is detected in the poorly-understood close circumstellar en- 
vironment (e.g. Jennings & Sada 1998, Verhoelst et al. 2006, 
Perrin et al. 2007). As we are limited by the transmission of 
Earth's atmosphere, our spectra do not cover the water bands 
which were discussed most (those with bandheads around 1.4 
and 1.9yum.) 

Regarding the identification of the feature observed here, 
H 2 is not a candidate since this molecule does not exhibit a 
clear band at this particular wavelength range. CN provides an 
excellent match to the bandhead, and by a careful choice of 
temperature (2 500K) and column density (6xl0 18 cm -2 ), it is 
possible to model this absorption feature with a single slab of 
CN, without the introduction of unobserved features at the other 
wavelengths of CN opacity (see the upper panel of Fig.[9]for the 
model fit). 

Unfortunately, the same basic modelling approach does not 
work for the emission spectra: the prediction of the bluest fea- 
ture is, in the optically thin case, far too narrow (lower panel of 
Fig. [9]) and any combination of temperature and column density 
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Fig. 9. Upper panel: The photospheric model spectrum (black), 
the model + CN slab spectrum (green), the predicted inband 
intensities (blue) and the observations (red). Lower panel: The 
emission spectrum of optically thin CN at 2 500 K in green, the 
observations in red at a radius of 6 R* and the predicted inband 
intensities in blue. Both plots use the same temperature and col- 
umn density for the CN slab. 



that matches this broad feature, predicts a strong emission fea- 
ture in the K band starting around 1.9 /mi, but is not observed. 
Neither dumpiness nor a difference between excitation temper- 
ature and source-function temperature or a source of continuous 
background emission help to solve this problem. Obscuration of 
the feature at 1 .9 /mi is in theory possible with cold water vapour 
(cfr. the discussion on the identification of the observed absorp- 
tion features at those wavelengths), but that would require a large 
water column density and a very particular geometry. 

An important aspect of our spectro-imaging observations is 
that the angular resolution decreases by a factor two in the K 
band compared to the / band (where the plume is most visi- 
ble). In the L-R deconvolution process, the angular size of the 
bright photosphere is recovered properly (Fig. [6]), but the de- 
tection of the plume is made more difficult by its much lower 
surface brightness (typically several hundred times fainter than 
the photosphere). We therefore cannot exclude that our marginal 



detection of the plume in the K band is caused by a lack of re- 
solving power, that results in a lower sensitivity of the L-R de- 
convolved images to extended sources. As a consequence, as we 
do not observe the expected emission feature around 1.9 /mi, we 
cannot confirm our idenfication of CN, but cannot propose a vi- 
able alternative either. 



4.2. Previous high angular resolution observations 

Thanks to its large angular diameter and correlatively high 
brightness, Betelgeuse has been an excellent target for high res- 
olution observations using a broad palette of techniques: di- 
rect full-pupil imaging, aperture masking, speckle interferome- 
try and long-baseline interferometry. This short summary of the 
existing observational reports is in no way complete, but is in- 
tended to give an overview of the detections of asymmetries and 
surface features on Betelgeuse. Such asymmetries could point at 
a link between the photosphere and the envelope of the star as 
observed in our NACO images. 

Early observations in the visible using the rotation shear- 
ing interferometer (Roddier & Roddier 119831 [T985; Roddier et 
al. 119861 ) revealed that in November 1980 and February 1982, 
asymmetric emission was present around Betelgeuse, and they 
proposed that it was caused by scattering in a dusty layer close 
to the star. The presence of two companions was suggested by 
Karovska et al. ( 1986 ), but our NACO observations do not show 
any point- like source down to a magnitude difference > 5 in the 
JHK band over the field of view of our observations (0.4 x OA"). 
Observations in the visible obtained at the same epoch (February 
1981) with the Mayall telescope by Christou et al. (1988) show 
the presence of PSF- subtraction residuals in the southwestern 
quadrant with a typical radius of the order of 100 mas, a size and 
position comparable to that of the plume detected with NACO. 
Hebden et al. (11987b observed Betelgeuse in the Ha line with 
the MMT in speckle interferometry at two epochs: December 
1983 and November 1985. They identify mostly no asymmetry 
at both epochs up to a radius of « 95 mas, while slightly larger 
extensions appear in the east- west and north- south directions at 
larger radii. In February 1989, Buscher et al. (119901) observed 
Betelgeuse using the WHT equipped with a non-redundant aper- 
ture mask and three narrow-band filters in the visible around 
X = 0.7 /mi. These authors detected an asymmetry on the sur- 
face of the star, that they interpret as a hot spot located at a PA 
of « 280°, with a relative flux contribution of 10-15%. Using the 
same instrument and technique, Wilson et al. (1992 ) re-observed 
Betelgeuse in January 1991, and also detected asymmetries on 
the disk of the star and a northeast- southwest elongation differ- 
ent from to the 1989 epoch. 

Gilliland & Dupree (1996) resolved the chromosphere of 
Betelgeuse at ultraviolet wavelengths using the Hubble Space 
Telescope equipped with the Faint Object Camera (FOC), 
and also obtained resolved spectra with the Goddard High 
Resolution Spectrograph (GHRS). These data were obtained in 
March 1995 and clearly show the resolved chromosphere in UV 
continuum, extending to more than twice the radius of the visible 
photosphere. The deconvolved FOC images consistently show 
an unresolved bright region located in the southwestern quadrant 
of the disk, where we observed the plume with NACO. Gilliland 
& Dupree ( 1996) estimate that the spot temperature is > 200 K 
hotter than the surrounding chromosphere. The chromosphere 
observed at X - 255 nm with HST/FOC is on average much hot- 
ter (^ 5000 K) than the overall effective temperature of the star 
(« 3600 K, Perrin et al. 120041) . 
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In October 1995, Burns et al. (fT997b used the COAST in- 
terferometer at A - 0.80 yum to observe Betelgeuse and found 
no asymmetry of the stellar disk itself, but their reconstructed 
image (their Fig. 4) shows an extension along the northeast- 
southwest direction (PA « 200°), aligned with the NACO plume. 
Further observations obtained in November 1997 with COAST 
and the WHT telescope at visible and near infrared wavelengths 
(0.70, 0.91 and 1.29 jim) by Young et al. (f2000b also showed an 
asymmetry of the light distribution on the disk, that tends to de- 
crease as wavelength increases. Interestingly, their 1.29 /mi ob- 
servations do not show any asymmetry, while our observations 
with NACO at 1.28 /mi still show the presence of the southwest- 
ern plume. Tuthill et al. (119971) also used the WHT telescope in 
aperture masking mode around A = 0.70 yum between January 
1992 and December 1993. Depending on the epoch, they identi- 
fied two to three hot spots with changing positions, and relative 
flux contributions of up to 20%, on the surface of Betelgeuse. 

Betelgeuse was observed in early 1998 by Hinz et al. ([1998 ) 
using nulling interferometry in the mid-infrared (A = lOyum) 
with two mirrors of the Multiple Mirror Telescope separated by 
a baseline of 5 meters. This interferometric technique allows to 
cancel out the light from the unresolved star, leaving as a resid- 
ual the emission from the circumstellar environment. Their re- 
sulting nulled image (their Fig. 3) shows a very extended (« 5" 
in diameter) asymmetric dusty envelope, with extensions in the 
southwestern, northwestern and northeastern directions, similar 
to the positions of the three plumes we observed with NACO. 
Smith et al. (12009b recently found a mildly clumpy, spherical 
shell around Betelgeuse in CO emission, from spatially resolved 
spectroscopy at A « 4.6/mi. Some interesting structures were 
detected « 1" west of the star, but their identification to a mass 
loss event or envelope asymmetry is unclear. 

The single-mode two-telescope interferometric observations 
obtained by Perrin et al. (12004b and Chagnon et al. (12002b re- 
spectively in the K and L bands provided high accuracy pho- 
tospheric angular diameters of Betelgeuse (#ldk = 43.7 ± 
0.10 mas, #udl = 42.17 ± 0.05 mas). However, due to the ab- 
sence of closure phase measurements, it was not possible to 
check for the presence of asymmetry on the stellar disk. New 
observations of Betelgeuse with the three-telescope beam com- 
biner IONIC installed at the IOTA interferometer were obtained 
in October 2005 and March 2006 by Haubois et al. (12006b in 
the infrared H band. They point at the presence of one very 
faint hot spot (relative flux 0.5%) in the northwestern quadrant 
of Betelgeuse (through model-fitting), or alternatively at two 
hot spots close to the center of the disk (through image recon- 
struction). Spectrally dispersed interferometric observations of 
Betelgeuse in the K band (obtained in January 2008) are reported 
by Ohnaka et al. d2009b . These authors find a centrally symmet- 
ric emission on the apparent disk of the star in the continuum, but 
an inhomogeneous velocity field in the 2.303 yum CO line, with 
a North-South preferential direction. This direction appears con- 
sistent with the preferential axis of the main plume we detected 
with NACO. 

At thermal infrared wavelengths (A = 11.15yum), Tatebe et 
al. d2007b identified an asymmetry on the disk of Betelgeuse 
from three-telescope interferometric observations with the ISI 
instrument obtained in November and December 2006. They in- 
terpreted this asymmetry as the signature of a hot spot located on 
the southern limb of the star at a PA of 170 - 180°. This position 
is compatible with the PA of the plume we detect with NACO. 
Perrin et al. d2007b observed Betelgeuse with the two-telescope 
instrument VLTI/MIDI in the N band (A = 7.5 - 13.5 /mi), de- 



tecting a molecular layer (the MOLsphere) containing H2O, SiO 
and AI2O3 at a radius of « 1.4 R* from the star. 

4.3. Rotation 

The rotation of a supergiant star as Betelgeuse is expected to be 
very slow, but it can have non-negligible effects on the stellar at- 
mosphere structure. Uitenbroek et al. (1998 ) presented a detailed 
analysis of the HST/GHRS spectra first reported by Gilliland 
& Dupree ( 1996). They showed that chromospheric emission in 
the UV emission lines of Mg II extends to « 270 mas from the 
center of the star, twice further than the UV continuum. In ad- 
dition, thanks to a scanning of the instrument over the stellar 
disk, they could determine the projected rotational velocity of 
(vsin/ = 5km/s), the PA of its polar axis (PA « 235°) and the 
inclination of the polar axis on the line of sight of (i = 20°). 
From spectroscopy and theoretical considerations, Gray ( 2000) 
obtained a similar rotational velocity. This position of the polar 
cap is compatible with the hot spot visible in the FOC images 
and the PA of the southwestern plume observed in our NACO 
images (PA= 200 ± 20°) 

4.4. Origin of the southwestern plume 

The presence of the bright southwestern plume detected with 
NACO implies that the spherical symmetry observed on the in- 
frared photosphere of Betelgeuse is not preserved in its close 
environment. Among the plausible mechanisms to explain its 
presence, two appear as promising: convection and rotation. 
Following the proposal by Schwarzschild dl975b and the sim- 
ulation results by Freytag & Hofner (2008), a first hypothesis 
is that the very large convection cells enhance locally the mass- 
loss from the star, creating the observed plumes of molecules. 
However, although many studies have shown that the photo- 
sphere of Betelgeuse presents a significant degree of asymme- 
try, variable with time, the link with the chromosphere and the 
envelope remains to be established. A second hypothesis is that 
mass loss may be enhanced above the polar cap of the star, cre- 
ating the southwestern plume of Betelgeuse. This possibility is 
supported by the HST imaging and spectroscopy by Uitenbroek 
et al. (11998b , as it reveals the presence of a hot spot that could 
be located at the pole of the star. An interaction between this hot 
spot and the extended chromosphere of the star could result in 
a preferential production of molecules above the pole and the 
creation of a plume. 



5. Conclusion 

From diffraction-limited observations in ten narrow-band filters 
spread over the JHK spectral domain, we find evidence for the 
presence of a complex CSE around Betelgeuse. The most con- 
spicuous feature of this envelope is a bright plume extending in 
the southwestern quadrant up to a minimum distance of six pho- 
tospheric radii. Its SED in emission appears complementary to 
a strong absorption feature present in the / band ( 1.08-1. 24 //m) 
on the photospheric spectrum. This points at the presence of a 
molecular layer, possibly CN, although the formal identification 
of the spectral signature of this molecule in the K band is not 
established. Two hypotheses could explain the formation of the 
plume: convection causing an enhanced mass loss above a large 
upwards moving convective cell, or rotation, through the pres- 
ence of a hot spot at the location of the polar cap of the star. 
Further spectro-imaging observations resolving the photosphere 
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of Betelgeuse (e.g. by interferometry) should allow to discrimi- 
nate between these possibilities, if a link can be established be- 
tween photospheric hot spots (characteristic of convection) and 
the detected plume. 
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Fig. 2. Average NACO images of Betelgeuse (left image of each column) and the PSF calibrator Aldebaran (right image of each 
column) in 10 narrow-band filters from 1.04 to 2.17/zm. The color scale is a function of the square root of the flux, normalized to 
the maximum and minimum value of each image, in order to emphasize the fainter parts of the images. The white disks in the lower 
right corners of the 1.04//m images represent the photospheric angular diameters of Betelgeuse and Aldebaran. 
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Fig. 6. Lucy-Richardson deconvolved images of Betelgeuse (30 iterations). For each wavelength, the left image is the full-dynamic 
linear color scale image normalized to the maximum and minimum values. The 10 contours plotted in red correspond to flux 
levels of 1 000, 2000, 4000, 8 000, 16000, 32000, 64000, 128000 and 256000 W.m~ 2 .fim' 1 .sr' 1 . The right image uses a square 
root color scale with narrower color cuts (including 95% of the pixels) to emphasize the faint extensions of the CSE. The green 
circle corresponds to the 43.7 mas photospheric angular diameter of the star as measured by Perrin et al. (2004) from K band 
interferometric observations. 
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05 


597 


Aldebaran 


NB1 


24 


• 


64 


7 


2 


15000 


0.73 


1 


.33 


046 


2009- 


-01 


-04T02 


14 


22 


059 


Aldebaran 


NB1 


09 


• 


64 


7 


2 


15000 


0.64 


1 


.33 


047 


2009- 


-01 


-04T02 


16 


38 


434 


Aldebaran 


NB1 


08 


• 


64 


7 


2 


15000 


0.57 


1 


.33 


048 
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-01 


-04T02 


18 


55 
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04 


• 


64 


7 


2 
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0.74 


1 


.33 


049 
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22 


17 
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12 


• 


64 


7 


2 


15000 


0.69 


1 


.33 


050 


2009- 


-01 


-04T02 


24 


33 


229 


Aldebaran 


NB2 


17 


• 


64 


7 


2 


15000 


0.60 


1 


.33 


051 


2009- 


-01 


-04T02 


32 


50 


331 


Betelgeuse 


NB1 


64 


• 


64 


7 


2 


15000 


0.57 


1 


.24 


052 


2009- 


-01 


-04T02 


34 


58 


410 


Betelgeuse 


NB1 


75 


• 


64 


7 


2 


15000 


0.59 


1 


.24 


053 


2009- 


-01 


-04T02 


37 


17 


272 


Betelgeuse 


NB1 


28 




64 


7 


2 


15000 


0.60 


1 


.23 


054 


2009- 


-01 


-04T02 


39 


29 


761 


Betelgeuse 


NB1 


26 




64 


7 


2 


15000 


0.65 


1 


.23 


055 


2009- 


-01 


-04T02 


41 


46 


478 


Betelgeuse 


NB1 


24 




64 


7 


2 


15000 


0.63 


1 


.22 


056 


2009- 


-01 


-04T02 


44 


04 


635 


Betelgeuse 


NB1 


09 




64 


7 


2 


15000 


0.65 


1 


.22 


057 


2009- 


-01 


-04T02 


46 


21 


058 


Betelgeuse 


NB1 


08 




64 


7 


2 


15000 


0.75 


1 


.22 


058 


2009- 


-01 


-04T02 


48 


37 


938 


Betelgeuse 


NB1 


04 




64 


7 


2 


15000 


0.65 


1 


.21 


059 
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-01 


-04T02 


51 


02 
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12 




64 


7 


2 


15000 


0.59 


1 


.21 
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Table 1. continued. 



050 


2009-01-04T02:53: 19.617 


Betelgeuse 


NB2.17 


• 


64 


7.2 


15000 


0.61 


1.21 


051 


2009-01-04T03:ll:04.218 


310ri 


NB1.64 


o 


64 


7.2 


15000 


0.69 


1.09 


052 


2009-0 1-04T03: 13:24.544 


310ri 


NB1.26 


o 


64 


7.2 


15000 


0.69 


1.09 


053 


2009-0 1-04T03: 15:41.604 


310ri 


NB1.08 


o 


64 


7.2 


15000 


0.73 


1.09 


054 


2009-01-04T03:18:05.278 


310ri 


NB2.17 


o 


64 


7.2 


15000 


0.79 


1.09 



a The "•" symbol indicates that the neutral density filter was used ("o" indicates it was not used). 
b DIT is the detector integration time. 

c cr is the DIMM seeing at Paranal in the visible and AM is the airmass. 



